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Abstract 

MUC5AC is a well-known gastric differentiation marker, which has been frequently used for the classification of stomach 
cancer. Immunohistochemistry revealed that expression of MUC5AC decreases accompanied with increased malignant 
property of gastric mucosa, which further suggests the importance of MUC5AC gene regulation. Alignment of the 5'- 
flanking regions of MUC5AC gene of 13 mammal species denoted high homology within 200 bp upstream of the coding 
region. Luciferase activities of the deletion constructs containing upstream 451 bp or shorter fragments demonstrated that 
15 bp region between —111 and —125 bp plays a critical role on MUC5AC promoter activity in gastrointestinal cells. We 
found a putative Gli-binding site in this 15 bp sequence, and named this region a highly conserved region containing a Gli- 
binding site (HCR-Gli). Overexpression of Gli homologs (Glil, Gli2, and Gli3) clearly enhanced MUC5AC promoter activity. 
Exogenous modulation of Glil and Gli2 also affected the endogenous MUC5AC gene expression in gastrointestinal cells. 
Chromatin immunoprecipitation assays demonstrated that Glil directly binds to HCR-Gli: Gli regulates ML/C5;AC transcription 
via direct protein-DNA interaction. Conversely, in the 30 human cancer cell lines and various normal tissues, expression 
patterns of IVIUCSAC and Gli did not coincide wholly: MUC5AC showed cell line-specific or tissue-specific expression 
whereas Gli mostly revealed ubiquitous expression. Luciferase promoter assays suggested that the far distal MUC5AC 
promoter region containing upstream 4010 bp seems to have several enhancer elements for gene transcription. In addition, 
treatments with DNA demethylation reagent and/or histone deacetylase inhibitor induced MUC5AC expression in several 
cell lines that were deficient in I\/1UC5AC expression. These results indicated that Gli is necessary but not sufficient for 
MUC5AC expression: namely, the multiple regulatory mechanisms should work in the distal promoter region of MUCSAC 
gene. 
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Introduction 

Although the mortality and incidence of gastric cancer has 
gradually fallen in the last several decades, it is still the fourth most 
common malignancy and the second leading cause of cancer- 
related death worldwide [1]. Gastric cancer mainly occurs from 
gastric mucosa with intestinal metaplasia, which is mostly caused 
by chronic infection of Helicobacter pylori [2,3]. Intestinal 
metaplasia in stomach is typically classified into two types: mixed 



gastric-and-intestinal type (incomplete type) and solely intestinal 
type (complete type) [4] . In the process of metaplastic change of 
gastric mucosa, induction of intestinal differentiation and loss of 
gastric differentiation occur simultaneously or asynchronously 

[5,6]. 

For induction of intestinal differentiation, we and other groups 
have been reported that Cdx must be an indispensable key 
molecule by regulating transcription of many intestinal marker 
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genes [7,8,9,10,11,12,13]. On the contrary, the molecular 
mechanism of loss of gastric genes in the process of metaplastic 
changi; in stomach remains poorly understood. It is probably due 
to inadequate identification of gastric marker genes, and also dtie 
to unsolved expression regulation of these gastric genes. Based on 
these backgrounds, we focused on MUC5AC, which is a well- 
established gastric marker gene [6,14] and is often used for the 
clinical assessment of gastric cancer [15]. There have been many 
previous reports investigating the expression of MUC5AC and 
prognosis of gastric cancer, but the association between expression 
of MIJC5AC and malignant potential of gastric cancer is still 
controversial [16,17,18,19]. Nevertheless, MUC5AC is one the 
most evident gastric marker clearly decreased in the process of 
intestinal metaplasia [14,20]. We believe elucidating the mecha- 
nism o{MUC5AC gene expression in gastrointestinal cells must be 
useful to understand the loss of gastric differentiation during the 
development of pre-malignant atrophic gastritis. 

Our aim of this study is to find a critical mechanism of 
MUC5AC expression regulation in human gastrointestinal cells. 
To date, HIF-la [21], Smad4 [21,22], Spl [22], GATA-4/-6, and 
HNF-1/-4 [23] are reported to activate murine MUC5AC. Spl 
and Gli are reported to enhance MUC5AC expression in human 
lung-epithelial and pancreatic cancer cells, respectively [24,25]. 
Contrastively, ATBFl is negatively regulate MUC5AC expression 
in human gastric cancer cells [26], In spite of these results, 
molecular mechanisms of MUC5AC expression in human 
gastrointestinal cells have not been fuUy elucidated. In the present 
study, we identified 15 bp sequence in human MUC5AC 5'- 
flanking region which plays an important role in MUC5AC 
expression in gastrointestinal cancer cells, and also found that the 
region contains a putative Gli-binding site which does not coincide 
with previous report [25]. Gli is one of transcriptional factors 
which has a DNA binding zinc finger domain [2 7] . In this study, 
we examined molecular roles of Gli on MUC5AC promoter in 
gastrointestinal cells. Our results demonstrated that MUC5AC 
expression depends on cooperative regulatory mechanism of Gli, 
some epigenetic modulation, and other factors in gastrointestinal 
cells. 

Materials and Methods 

Cell Culture 

Twenty gastric cancer cell hnes, ten colorectal cancer cell Unes, 
and two non-gastrointestinal cancer cell lines were maintained in 
high-glucose DMEM with 10% fetal calf serum (Gibco/Invitro- 
gen, Carlsbad, CA) at 37°C in a humidified 5% CO2 atmosphere. 
Names of used cell lines and histological types of gastric cancer cell 
hnes were described in our previous reports [6]. Human T98G 
and A 172 cell hnes were purchased from the RIKEN Bio 
Resource Center (Tsukuba, Japan). For the treatment with DNA 
demethylation reagent or histone deacetylase (HDAC) inhibitor, 5- 
Aza-2'-deoxycytidine (5-Aza-dC, Sigma-Aldrich) at 2 |J.g/ml and/ 
or trichostatin A (TSA, Sigma-Aldrich) at 25-1000 ng/ml were 
added to the culture medium. 

Tumor Samples 

For the advanced gastric cancer specimens, we randomly 
selected 89 gastric adenocarcinoma samples surgically resected at 
the Fujita Health University Hospital. For the early stage gastric 
cancer endoscopically resected, we selected 78 specimens banked 
at the University of Tokyo Hospital. This study was approved by 
the ethic committees of the University of Tokyo, and also by the 
institutional ethical review board for human investigation at Fujita 
Health University. According to the Declaration of Helsinki, 



written informed consents were obtained from all the study 
participants for use of resected sample in research. 

Immunohistochemistry 

Deparaflinization and endogenous peroxidase inactivation of 
cUnical tissues were performed as described previously [28]. For 
MUC5AC, hydrated heating in 1 mM EDTA buffer (pH 8.0) at 
120°C was then performed in a pressure cooker (Delicio 6L; T- 
FAL, RumUy, France) for 10 min for antigen retrieval. The 
primary immunostaining with anti-MUC5AC antibody (NCL- 
MUC-5AC, Novocastra, Newcastle-upon-Tyne, UK) at a 1:200 
dilution was applied for 12 hours at room temperature. After 
washing in PBS three times, the secondary immunostaining with 
Histofine Simple Stain MAX-PO(G) (Nichirei, Tokyo, Japan) was 
applied for 30 min at room temperature. Based on the evaluation 
by the two independent pathologists, ratios of MJ7C5/lC-positive 
cells were classified into four classes: 1) <10%, 2) S10% and < 
50%, 3) >50% and <80%, and 4) >80%. 

For Glil immunostaining, hydrated heating in 10 mM citrate 
buffer (pH 6.0) was done in a pressure cooker for 10 min for 
antigen retrieval. The sections were then incubated for 1 h at 
room temperature with antihuman Glil antibody (H-300, sc- 
20867, Santa Cruz Biotechnology) at a 1:500 dilution. For the 
ampKfication of signals, anti-rabbit antibody (Dako, Hamburg, 
Germany) was apphed to the slides for 30 min at room 
temperature. This was followed by incubation with FITC- 
conjugated phenol (fluorescyl-tyramide; Dako) for 1 5 min at room 
temperature and incubation with anti-FITC antibody conjugated 
to HRP (Dako) for 15 min at room temperature was then done. 

Plasmld constructions 

The primer sequences used in plasmid constructions are listed in 
Table 1 . To construct the vectors for luciferase reporter assays, the 
2000 bp and 1433 bp of MUC5AC promoter region were 
amplified from TIG-112 genome using the primers MUC5A- 
Cup-F-01/MUC5ACup-R-01 and MUC5ACup-F-02/MUC5A- 
Cup-R-01. The amplified products were cloned into pT7blue-T- 
vector (Novagen, Darmstadt, Germany), and then digested with 
Sail and BamHI. The excised 2 kb and 1.4 kb DNA fragments 
were inserted into the Xhol/Bglll site of pGL4. 1 2 (Promega, 
Madison, WI, USA) to generate pGL4.12-hMUC5ACup2000bp 
and pGL4.12-hMUC5ACupl433bp. Ncol-Hindlll fragment 
from pGL4.12-hMUC5ACup2000bp was inserted into the 
EcoRV/Hindlll site of pGL4.12 to generate pGL4.12-hMUC5A- 
Cup451bp. The 280 bp, 190 bp, 158 bp, 150 bp, 138 bp, 125 bp 
and 110 bp oi MUC5AC promoter region were amplified from 
pGL4.12-hMUC5ACup451bp as a template using primers 
MUC5ACup-280-XhoI, MUC5ACup-190-XhoI, MUC5ACup- 
158-XhoI, MUC5ACup-150-XhoI, MUC5ACup-138-XhoI, 
MUC5ACup-125-XhoI, MUC5ACup-l 10-XhoI and pGL412- 
Hindlll, digested with Hindlll and Xhol, and inserted into the 
Hindlll/Xhol site of pGL4.12 to generate pGL4.12-hMUC5A- 
Cup280bp, pGL4.12-hMUC5ACupl90bp, pGL4.12-hMUC5A- 
Cupl58bp, pGL4.12-hMUC5ACupl50bp, pGL4.12-hMUC5A- 
Cupl38bp, pGL4.12-hMUC5ACupl25bp and pGL4.12- 
hMUC5ACupll0bp. The 4010 bp to 1432 bp upstream of 
MUC5AC promoter region was amplified from genome of TIG- 
1 12 genome using the primers MUC5ACup-SacI-F and MUC5A- 
Cup-PciTR, digested with Sacl and Pcil, and inserted into the 
Sacl/Pcil site of pGL4.12-hMUC5ACup2000bp to generate 
pGL4.12-hMUC5ACup4010bp. The 3000 bp to 1432 bp up- 
stream of MUC5AC promoter region was amplified from TIG- 
112 genome using the primers MUC5ACup-3000-NheI-F and 
MUC5ACup-PciI-R, digested with Nhel and PcU, and inserted 
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into the Nhel/Pcil site of pGL4.12-liMUC5ACup2000bp to 
generate pGL4. 12-hMUC5ACup3000bp. To delete tlie 15 bp in 
the MUC3AC promoter region (125 bp to 1 11 bp upstream), 
PCR-based mutagenesis was performed using primers MUC5A- 
Cup-15delete-F and MUC5ACup-15delete-R and cloned into 
pGL4.12 to generate pGL4.12-hMUC5ACup451bp-deltal5bp. 
The plasmid was confirmed by sequence analysis. The 284 bp 
Sbn-Hindlll region of pGL4. 12-hMUC5ACup4010bp was re- 
placed by the 269 bp Sbll-Hindlll fragment of pGL4.12- 
hMUC5ACup451 bp-delta 15bp to generate pGL4. 12-hMUC5A- 
Cup40 1 Obp-delta 1 5bp. 

Glil and Gli3 cDNAs were purchased from Open Biosystems 
(clone ID 3531657 and 40125719, respectively, Huntsville, AL,) 
and Gli2 cDNA (pENTR223.1-hGli2) was purchased from 
DNAFORM K.K. (clone ID 100069128, Yokohama, Japan). 
EcoRTXhoI fragment containing Glil cDNA was inserted into 
the EcoRI/XhoI site of pMXs-IP to generate pMXs-Glil-IP. 
EcoRTEcoRV fragment containing Gli3 cDNA was inserted into 
the EcoRI/PacI site of pMXs-IP to generate pMXs-Gli3-IP. 
Xbal-NotI fragments from pMXs-Glil-IP and pMXs-Gli3-IP 
were inserted into the Nhel/NotI site of pcDNA3.1(4-) to generate 
pcDNA3.1(+)-Glil and pcDNA3.1(+)-Gli3. BsrGI-Avrll fragment 
from pENTR223.1-hGli2 was inserted into the Acc65I/XbaI site 
of pcDNA3.1(+) to generate pcDNA3.1(+)-Gli2. 

Luciferase reporter assay 

Cells were cultured on 96-well plates and transiently transfected 
with mixtures of Renrlla luciferase control vector pGL4.74 
(Promega) (15 ng) and the firefly luciferase reporter plasmids 
(150 ng) by using Lipofectamine and Lipofectamine PLUS 
(Invitrogen). To examine the effect of Gli on the promoter 
activities, pcDNA3.1(4-)-Glil, -Gli2, and -Gh3 plasmids were 
transfected with pGL-based plasmids at the same time. Luciferase 
assays were performed at 24 or 48 h post-transfected using the 
Dual luciferase reporter assay system (Promega). Luciferase 
activities were measured by an AutoLumat Plus LB953 (Berthold, 
Bad WUdbad, Germany). Luciferase activity was normalized to a 
ReniUa control and the results are shown as mean ± SD from at 
least three independent experiments. 

Production and infection of retrovirus vectors 

pMXs-Glil-IRES-Puro, pMXs-Gli2-IRES-Puro and pMXs- 
Gli3-IRES-Puro and control vector (pMXs-IRES-Puro) were 
transfected with pCAG-VSVG into the PLAT-GP packaging cell 
line for retrovirus vector production (Cell Biolabs, Inc. San Diego, 
CA) using Lipofectamine and Lipofectamine PLUS (Invitrogen). 
Virus-containing supernatants were harvested at 24 h, 48 h and 
72 h post-transfection, pooled, and filtered though a 0.45-|J,m filter 
(MrUipore, Bedford, MA). To assess the effects of overproduction 
of Gli, cells were infected by adding retrovirus-containing 
supernatant at a final concentration of 8 \ig/ ml Polybrene (Sigma) 
and were selected over 3 days with puromycin (Sigma). 

Transfection 

Cells were transiently transfected with pcDNA3. l(4-)-Glil, 
pcDNA3.1(+)-Gli2, pcDNA3.1(+)-Gli3 and pcDNA3.1(+) plasmids 
for 24-48 h using Lipofectamine and Lipofectamine PLUS, and 
were used for RT-PCR and luciferase reporter assay. 

RT-PCR 

Total cellular RNA was prepared using the Isogen RNA 
isolation reagent (Wako Pure Chemical Industries, Osaka, Japan) 
as previously reported [29]. RT-PCR was performed via a 



Superscript One-Step reaction using the Platinum Taq (Invitro- 
gen). The primer sequences used for RT-PCR are shown in 
Table 2. RNA was reverse-transcribed for 30 min at 50°C, and 
after an initial denaturation at 94°C for 3 min, cDNA amplifica- 
tion procedures were performed as follows: for MUC5AC, 35 
cycles of 94°C for 30 sec, 64°C for 1 min, and 72°C for 1 min; for 
GAPDH, 30 cycles of 94°C for 30 sec, 60°C for 1 min, and 72°C 
for 1 min; for Glil and Gli2, 35 cycles of 94°C for 30 sec, 60°C 
for 1 min, and 72°C for 1 min; for GU3 , 40 cycles of 94°C for 
30 sec. A commercial RNA panel. Human Total RNA Master 
Panel II, was purchased from Clontech Laboratories (Palo Alto, 
CA, USA). 

siRNA 

SH-IO-TC ceUs were transfected with siGENOME SMART- 
pool siRNA against human Glil and/or Gli2 (M-003896-00-0005 
and M-006468-02-0005, respectively; Dharmacon, Lafayetta, CO, 
USA) using Lipofectamine RNArMAX (Invitrogen) according to 
the manufacture's instructions. RNA was extracted from the cells 
48 h after siRNA transfection and used for RT-PCR. 

Chromatin immunoprecipitation (ChIP) assay 

Chip analyses were performed using a ChIP assay kit (Upstate 
Biotechnology Inc., Lake Placid, NY) according to the manufac- 
ture's instructions. For crosslinking, cells (2x10**) were incubated 
at 37°C for 15 min in PBS containing 1% formaldehyde. Cells 
were collected after two washes with ice-cold PBS and subjected to 
centrifugation for 5 min at 700 xg. Cells pellets were resuspended 
400 10.1 sodium dodecyl sulfate (SDS) lysis buffer (ChIP assay kit. 
Upstate) containing protease inhibitors (0. 1 M phenylmethylsulfo- 
nylfluoride and 4 |j.l protease inhibitor cocktail (Sigma)), and 
incubated for 10 min on ice. DNA was sonicated (setting 5, Handy 
Sonic, model UR-20P; Tomy Seiko, Co., Ltd., Tokyo, Japan) and 
then subjected to centrifugation at 13,000 rpm for 10 min. 
Supernatants were diluted 10-fold with Chip dilution buffer (ChIP 
assay kit. Upstate). 1% of the supernatant was retained as the 

Table 1. Primers used for plasmid constructions. 



Name 


Sequence 


MUC5ACup-F-01 


5 '-attcattcacccattcactcactcattcac-3 ' 


MUC5ACup-F-02 


5'-tgccacatgtgaagtgctctttctctaggc-3' 


MUC5ACup-R-01 


5'-tgtgtggacggcggggaagagtgccctgtc-3' 


MUC5ACup-280-Xhol 


5'-taggctcgagtatgtggggaggacccctgc-3' 


MUC5ACup-190-Xhol 


5 '-cgccctcgagcggccgctggccagcccgca-3 ' 


MUC5ACup-158-Xhol 


5'-gagcctcgagtgtttacttgggtgaggggg-3' 


MUC5ACup-110-Xhol 


5 '-gcccctcgaggtgaagcacggggctggagc-3 ' 


MUC5ACup-150-Xhol 


5'-actgctcgagtgggtgagggggaaccacag-3' 


MUC5ACup-138-Xhol 


5'-ggtgctcgagaaccacaggccccgccctgc-3' 


MUC5ACup-125-Xhol 


5'-cacactcgaggccctgcccacccacgtgaa-3' 


pGL412-Hindlll 


5 '-ccggattgccaagcttggcc-3 ' 


MUC5ACup-Sacl-F 


5'-tgcccacagagctcagaaacaaggc-3' 


MUC5ACup-Pcil-R 


5'-agagcacttcacatgtggcaggagt-3' 


MUC5ACup-3000-Nhel-F 


5'-acttgctagctcatacactcattcacttat-3' 


MUC5ACup-15delete-F 


5'- 

gagggggaaccacaggccccgtgaagcacggggctggagc- 
3' 


MUC5ACup-15delete-R 


5'-gctccagccccgtgcttcacggggcctgtggttccccctc-3' 


doi:l 0.1 371 /journal.pone.Ol 061 Oe.tOOl 
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Table 2. Primer pairs, annealing temperatures (Tm), and product sizes (Length) for the 5 genes analyzed by RT-PCR. 



Genes 




Forward (F) and reverse (R) primer sequences 


Tm(°C) 


Length 


MUC5AC 


F 


5 '-accggtgccacatgacggac-3 ' 


64 


396 




R 


5'-acgtggccgcctcacacgtg-3' 






Gli1 


F 


5 '-gaaggagttcgtgtgccact-3 ' 


60 


491 




R 


5'-gtctgctttcctccctgatg-3' 






Gli2 


F 


5 '-gcaacaaagccttctccaac-3 ' 


60 


430 




R 


5 '-atctccacgccactgtcatt-3 ' 






Gli3 


F 


5'-aaagcaaacaggagcctgaa-3' 


60 


401 




R 


5 '-ggaatgcgttctgttttggt-3 ' 






GAPDH 


F 


5 '-accacagtccatgccatcac-3 ' 


60 


423 




R 


5'-tccaccaccctgttgctgta-3' 







doi:l 0.1 371/journal.pone.Ol 061 06.t002 



input, and the rest was then subjected to immunoprecipitations. 
Immunoprecipitations were performed overnight at 4°C with 3 |J,g 
of anti-GHl antibody (H-300, sc-20687; Santa Cruz Biotechnol- 
ogy, Santa Cruz, CA) and nonimmunized rabbit IgG whole 
molecule (sc-2027; Santa Cruz Biotechnology). After reverse 
crosslinking, the obtained DNA was purified using PCR product 
purification kit (Qiagen, Valencia, CA, USA). Immunoprecipitat- 
ed DNA was analyzed by PCR using primers 5'-ctcggaaactgggctc- 
tacccgg-3' and 5'-gagctttttgtagccccagagctgg-3' to amplify a 
fragment of the MUC5AC promoter region. Primer pairs for the 
villin 1 promoter region were used as a negative control [8] . 

Results 

Decrease of MUCSAC expression accompanied with 
progression of gastric canceration 

To evaluate the expression of MUCSAC in gastric cancer, 
immunostaining was performed using the clinical specimen 
representing three-grade gastric epithelial cells from the view of 
tumorigenesis: non-tumorous but precancerous cells of atrophic 
mucosa around early gastric cancer (Fig. lA), malignant cells of 
early gastric cancer endoscopicaUy resected (Fig. IB), and 
malignant cells of advanced gastric cancer surgically resected 
(Fig. IC). Our results revealed that the stronger the malignant 
property of gastric epithelial cells is, the weaker the expression of 
MUCSAC is. In other word, expression of MUCSAC decreases 
accompanied with increased malignant property of gastric 
mucosa. 

Some previous studies including ours reported that expression of 
MUCSAC often attenuates during the development of pre- 
malignant intestinal metaplasia [6,14]. Our present data further 
showed that decrease of MUCSAC expression advances in the 
process of gastric canceration. From these results, we are 
convinced that it is important to elucidate the regulatory 
mechanism of MUCSAC gene expression for understanding 
tumorigenesis of gastric cancer. 

A highly conserved region containing a Gli-binding 
sequence (HCR-Gli) is present in the promoter of 
MUCSAC gene 

To examine the transcriptional regxilation of MUC5AC gene, 
the upstream S '-flanking regions oi MUC5AC gene of 13 mammal 
species retrieved from GenBank (dated 8 May, 2013) were aligned 
using ClustalW (ver 2.1) sequence alignment program [30]. The 



major transcription initiation site of human MUCSAC gene has 
been mapped to 48 bp upstream of the ATG translation start site 
[31]. Alignment of the human sequence with other mammal 
sequences denoted high homology among orthologs within 
~200 bp upstream of the ATG start codon (Fig. 2). 

To test the function of the conserved promoter regions 
identified above, a series of luciferase reporter constructs 
containing 4S1 bp or shorter fragments of the human MUCSAC 
promoter region were generated. Promoter activities of these 
deletion constructs were measured by luciferase reporter assay in 
the three M[/C5/lC-expressing cell lines derived from human 
gastrointestinal cancer: SW480, SH-IO-TC, and KE-39 [6]. As 
shown in Fig. 3A, the promoter activity of the reporter construct 
containing upstream 1 1 0 bp was much lower than those 
containing upstream 4S1, 280, 190, and 158 bp. 

To examine the essential region for the luciferase activity 
between — 1S8 to —110 bp more precisely, reporter constructs 
containing ISO, 138, and 12S bp upstream of ATG start codon 
were generated and examined their activities (Fig. 3B). Although 
the reporter constructs containing upstream 1S8, ISO, 138, and 
12S bp showed high luciferase activities, the construct containing 
upstream 1 10 bp denoted very low activities equal to a negative 
control. 

Next we made a deletion mutant construct containing upstream 
4S1 bp but lacking the region between — 12S and — 1 1 1 bp, and 
examined the promoter activity (Fig. 3C). Deletion of the IS bp 
led to a marked decrease in the luciferase activity, suggesting that 
the IS bp sequence between — 111 to — 12S bp is critical for the 
MUCSAC promoter activity. Based on the thorough literature 
search [32,33,34,3S,36,37,38], we found a putative Gli-binding 
site in this short IS bp sequence (S'-GCCCTGCCCACCCAC-3' 
shown in Fig. 3D). We could not find any other transcription 
factors which bind to this IS bp sequence by searching transcrip- 
tion factor database. Consequentiy, we named this IS bp sequence 
as a highly conserved region containing a Gli-binding site (HCR- 
Gli). 

Exogenous modulation of Gli affects the endogenous 
MUCSAC gene expression in gastrointestinal cells 

We examined whether Glil can upregulate MUCSAC tran- 
scription activity in gastrointestinal cell lines. We constructed a 
retrovirus vectors carrying the Glil gene, and transduced it into 
four gastric (SH-IO-TC, MKNl, NCI-N87, AZ-S21), one 
colorectal (WiDr), and three other cancer cell lines (MDA- 
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^ Precancerous epithelial g Malignant cells of early Q Malignant cells of 
cells of stomach mucosa gastric cancer " advanced gastric 

with atrophic gastritis endoscopically resected cancer surgically resected 




Value of MUC5AC expression 

Figure 1. Expression of IVIUC5AC in non-tumorous epitKielial cells of atrophic mucosa around early gastric cancer endoscopically 
resected (A), malignant cells of early gastric cancer endoscopically resected (B), and malignant cells of advanced gastric cancer 
surgically resected (C). Based on the percentages of cells with immunoreactivity of MUC5AC, the 78 non-malignant or malignant specimens (A, B) 
and 89 malignant specimens (C) were classified into four categories: 1) <10%, 2) >10% and <50%, 3) £50% and <80%, and 4) £80%. 
doi:1 0.1 371/journal.pone.01 061 Oe.gOOl 



MB435, T98G, A172). In the resulting stable transductants, the 
expression levels of MUC5AC, GUI, and GAPDH were analyzed 
by RT-PCR (Fig. 4A). Although there were two exceptions of AZ- 
521 and WiDr, MUC5AC was found to be universally upregulated 
by Glil overexpression. 

We next compared the three Gli homologs {GUI, GU2, and 
GU3) identified in human. To examine the elfect of these 
transcription factors on MUC5AC expression, AGS and SW480 
cell lines transiently overexpressing Glil, Gli2 and Gli3 were 
analyzed by RT-PCR (Fig. 4B). Glil strongly upregulated 
MUC5AC expression, which was consistent with the result of 
stable transductants in Fig. 4A. Gli2 and Gli3 also increased the 
MUC3AC expression, but their effects were weaker than Glil. 

To address whether GKl, Gli2, and Gli3 activate the MUC5AC 
promoter through the 15 bp of HCR-Gli, we performed the 
luciferase reporter assay using the construct containing upstream 
451 bp of the human MUG 5 AC gene and the deletion mutant 
lacking the HCR-GK (Fig. 4C). In the Gfi-transfected SW480, SH- 
10-TC, and KE-39 cells, all the three Gli homologs enhanced 
activation of the reporter construct containing upstream 45 1 bp of 
MUC5AC, but not of the deletion mutant lacking the HCR-Gli. 
These results indicate that Glil, Gli2, and Gli3 can efficiently 
activate MUC5AC transcription through the HCR-Gli. From 
thes(; rc-sults, speculate that Glil, Gli2 and Gli3 may play 
redundant roles but the effect of Glil is stronger than that of others 
upon the upregulation of MUC5AC expression. 



We next examined the effect of Gli knockdown by siRNA 
transfection on MUG5AC expression. After SH-IO-TC cells were 
transfected with Glil and/ or Gli2 siRNA, the expression levels of 
MUC5AC, GUI, Gli2, and GAPDH were analyzed by RT-PCR. 
As shown in Fig. 4D, double knockdown of Glil and Gli2 most 
efficiently decreased MUG5AC expression, whereas the effects of 
either Glil or Gli2 single knockdown were considerably weaker. 
These results suggest that GUI and Gli2 have redundant functions 
for MUC3AC expression, which is consistent with our above- 
mentioned speculation. 

Gli interacts to HCR-Gli in the MUC5AC promoter region 

To determine whether Glil directly binds to the MUG 5 AG 
promoter, chromatin immunoprecipitation (ChIP) assays were 
performed using anti-Glil antibody. As a negative control, the 
human viUinl promoter sequence was used because vilUnl is a 
typical intestinal marker gene regulated by Cdx [8], and also 
because exogenous modulation of Glil expression level did not 
change endogenous viUinI expression (our unpublished observa- 
tion). As shown in Fig. 4E, the HCR-Gli-containing fragment of 
MUG5AG promoter was co-immunoprecipitated with Glil in 
SW48() and KE-39 cells, whereas the upstream sequence oivilUnl 
promoter was not. Our results demonstrate that Glil activates 
MUG5AG gene transcription through direct interaction between 
Glil and the MUC3AG promoter region. 

Taking all the results from luciferase assays with upstream 
deletion constructs, overexpression and knockdown of Gli, and 
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-190 

human CAGACACACCCATCTCCGCCTGCCACCGGCCGCTGGCCAGCCCGCAGTGA GCACCC 

chimpanzee CAGACACACCCATCTCCGCCTGCCACCGGCCGCTGGCCAGCCCGCAGTGA GCACCC 

bonobo CAGACACACCCATCTCCGCCTGCCACCGGCCGCTGGCCAGCCCGCAGTGA GCACCC 

olivebaboon CAGACACACCCATCTCTGCCTGCCACCGGCCGCGGGCCAGCCCGCAGTGA GCACCC 

rhesusmonkey CAGACACACCCATCTCTGCCTGCCACCGGCCGCGGGCCAGCCCGCAGTGA GCACCC 

killerwhale CAGACACACCCATCCCTGCCTGTCTTCAGCTGCCGGCCAGCCCGTAGCCA GTGCCC 

dolphin CAGACACACCCATCCCTGCCTGTCTTCAGCTGCCGGCCAGCCCGTAGCCA GTGCCC 

horse CACCTGGGAG-ATAAC-ACTTGTCTCCAGCCGCCAGCCAGCTTGGAGTGA GTACCC 

dog CAGCTATGCCCGTCTCTGCCAGTCTTGGGCTGCTGGCCAGCCCCGGGTGA GTGCCC 

cat CGGCCACGCCCATCTCTGCCCGTCTCCGGCTGCTGGCCAGCCCCGGGCGA GTGCCC 

rat CCGGCACACCCATCCCTGACTGCT GCTATTACCCAGCCTGAAGTGACTGAGTGCCC 

panda CCGGCACACCCATCCCTGACTGCT GCTATTACCCAGCCTGAAGTGACTGAGTGCCC 

hamster CAGACACACCCACCTCTGACTGCT GCCATTGACTGGCCTAAAGTGACTGAGTACCC 

* * * ** * ** * * * *** 

-158 -150 -138 -125 -111 

human AC'tGTTTACT'tGGGTGAGGGGGkACCACAGGCCCCGCCCTGCCCACCCACGTGAAGCAC- 

Chimpanzee ACTGTTTACTTGGGTGAGGGGGAACCACAGGCCCCGCCCTGCCCACCCACGTGAAGCAC- 

bonobo ACTGTTTACTTGGGTGAGGGGGAACCACAGGCCCCGCCCTGCCCACCCACGTGAAGCAC- 

olivebaboon ACTGTTTACTCGGGTGAGGGGGAACCACAGGCCCTGCCCCGCCCACCCATGTGAAGCTC- 

rhesusmonkey ACTGTTTACTCGGGTGAGGGGGAACCACAGGCCCTGCCCCGCCCACCCATGTGAAGCTC- 

killerwhale ACTGTTTATCTGGGTG-AGGGGAACCACAGGCCCTGCCCCACCCTCCCACGTGAGGCCCG 

dolphin ACTGTTTATCTGGGTG-AGGGGAACCACAGGCCCTGCCCCACCCTCCCACGTGAGGCCCG 

horse ACTGTTTACCTGGGTG-AGGGGAAC-GCAG-CCCTGCCCCGCCCACC-ACGTGAAGAGGC 

dog ACTGTTTACCTGGGTG-GGGGGAACGGCAGGCCCCGCCCCACCCACCCACT 

cat ACTGTTTACCTGGGTG-AGGGGAACGGCAGGCCC-GCCCCGCCCACCCACGTGGAGCCAC 

rat ACTGTTTACCTGGGTGAGGGGGAACCACAGGCC-TGCCCCACCCACACACGTGAAGCCAC 

panda ACTGTTTACCTGGGTGAGGGGGAACCACAGGCC-TGCCCCACCCACACACGTGAAGCCAC 

hamster ACTGTTTACCTGGGTGAGGGG-AATGATAGGCC-TGCCCCGCCCACCCAAGGAAAGCCAC 
******** ***** *** ** ** ** **** *** * * 

TATA box TSS 

human -GGGGCTGGAGCCAGCTCTGGGGCTACAAAAAGCTCCTGCCACCTTGGGTCCCTCpTCAG 

chimpanzee -GGGGCTGGAGCCAGCTCTGGGGCTACAAAAAGCTCCTGCCACCTTGGGTCCCTCCTCAG 

bonobo -GGGGCTGGAGCCAGCTCTGGGGCTACAAAAAGCTCCTGCCACCTTGGGTCCCTCCTCAG 

olivebaboon -GGGGCCGGAGCCAGCTCTGGGGCTACAAAAGGCTCCTGCAACCTTGGGTCCCTCCTCAG 

rhesusmonkey -GGGGCCGGAGCCAGCTCTGGGGCTACAAAAGGCTCCTGCAACCTTGGGTCCCTCCTCAG 

killerwhale TGGGGCTGGAGCCAGCTCTGGGACTACAAAAGGCTCCTGCAGAGGCAGAGCCCTCCTAAG 

dolphin TGGGGCTGGAGCCAGCTCTGGGACTACAAAAGGCTCCTGCAGACGCAGAGCCCTCCTAAG 

horse -TGGGCTGGGGCCAGCTCTGGGGCTACAAAAGGCTCCTGCCGACTCGGGTGCCTCCTCAG 

dog — GGGCCTGAGCCAGCTCTGAGGCTACAAAAGGCTCCTGCCCACTCAGGTCCTTCTTCAG 

cat GGGGACCGGAGCCAGCTCTGAGGCTACAAAAAGCTCCCGGCCACTAGGGCTCTTCCTCAG 

rat TGGGGCTGGAGCCAGCTCTTGGACTACAAAAGGCTGCTGGTGACTAGAGTCTCTCCTCAG 

panda TGGGGCTGGAGCCAGCTCTTGGACTACAAAAGGCTGCTGGTGACTAGAGTCTCTCCTCAG 

hamster TGGGGCTGGAGCCAGTTCTCAGACTACAAAAGGCAGCTGGCCACTTGAGTCCCTCCTCAG 
** * * ***** *** * ******** ** * * ** * ** 

-1 

human AGGCTGCTGAGGGACAGGGCACTCTTCCCCGCCGTCCACAcAaTG 

chimpanzee AGGCTGCTGAGGGACAGGGCACTCTTCCCCGCCGTCCACACGATG 

bonobo AGGCTGCCGAGGGACAGGGCACTCTTCCCCGCCGTCCACACGATG 

olivebaboon AGGCTGCCGAGGGACAGGGCACTCTTCCCCGCCATCCACACGATG 

rhesusmonkey AGGCTGCCGAGGGACAGAGCACTCTTCCCCGCCGTCCACACGATG 

killerwhale AGGCCGCC CTGGGGACTTCTCCCC-CCACACACC--ATG 

dolphin AGGCCGCC CTGGGGACTTCTCCCC-CCACACACC--ATG 

horse AGGCCACCC CTGGGGACTTCTCCCCTTCCTCCACACGATG 

dog ACGCTGCCC CTGGGGACCGCTCTCCTCCATCCGCAGGATG 

cat AGGCTGCCC CCGGGGACTTCTCCCCTCCCTCTGCAGGATG 

rat AGGCTGCCAG TGGGGACATTTCCCCATGCTCCATAGCATG 

panda AGGCTGCCAG TGGGGACATTTCCCCATGCTCCATAGCATG 

hamster AGCCTGCCGG TGGGGACGTTTCCCCATGCTGCATAGCATG 

* * * * * ** ** ** 



Figure 2. Comparison of the 5'-upstream sequences of MUC5AC gene of 13 mammal species. Alignment was carried out using the 
ClustalW (ver 2.1). The translational start site (TSS, arrow) and the translational starting codon (ATG, bold) are indicated. Asterisks (*) represent exact 
matches in all the sequences. 
doi:1 0.1 371/journal.pone.01 061 06.g002 
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Figure 3. Luciferase reporter anaiyses of a series of A7(/(r5/ICpromoter constructs in gastrointestinai ceii lines and identified tiighly 
conserved sequence similar to Itnown Gii-binding sequences. (A-C) Data represent the mean of luciferase activities of SW480, SH-10-TC, and 
KE-39 cells measured at 24 h after transfection. The error bars showed the standard deviation of the results from three independent experiments. (D) 
Alignment of the highly conserved sequence in MUC5AC promoter region with previous reported sequences of Gii-binding site. 
doi:1 0.1 371/journal.pone.01 061 06.g003 
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Figure 4. Effects of exogenous modulation of Gli on /Wi/C5/lCgene expression. (A) RT-PCR detecting MUC5AC mRNA in 5 gastrointestinal 
and 3 non-gastrointestinal cell lines infected with retroviral vector encoding GUI gene. (B) RT-PCR detecting MUC5AC mRNA in AGS and SW480 cells 
transfected with pcDNAS.IW (vec), pcDNA3.1(4-)-Gli1, pcDNAS.I (4-)-Gli2 or pcDNA3.1 (4-)-Gli3. (C) Luciferase reporter analysis of MUC5AC promoter 
constructs in SW480, SH-10-TC, and KE-39 cells transfected with pcDNA3.1W (vector), pcDNA3.1(4-)-Gli1, pcDNA3.1(-l-)-Gli2 or pcDNA3.1 (-l-)-Gli3. (D) 
MUC5AC expression analyzed by RT-PCR using SH-10-TC cells transfected with control siRNA, Glil sIRNA alone, GII2 siRNA alone, or Glil and GII2 
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siRNA. (E) Chip analyses of the human MUCSACqene in SW480 and KE-39 cells. Anti-Glil antibody (Glil-Ab) and non-immunized rabbit IgG whole 
molecule were used for immunoprecipitation. PCR was performed with primers recognizing the HCR-Gli-containing fragment MUC5AC promoter or 
recognizing the promoter region of villinl gene. 
doi:l 0.1 371/journal.pone.01 061 06.g004 



Chip analyses into consideration, we concluded that Gli enhances 
the MUC5AC gene transcription via direct protein-DNA interac- 
tion on the HCR-Gli region. 

Gli is not sufficient for transcription of MUC5AC gene in 
gastrointestinal cells 

To examine the relationship between endogenous expression 
levels Gli and MUC5AC in gastrointestinal cells, RT-PCR was 
performed using the twenty gastric, ten colorectal, and two non- 
gastrointestinal cancer cell lines (Fig. 5A). AU the cell lines with 
endogenous MUC5AC expression were found to be positive for 
GUI and/or Gli2 expression. Conversely, there were several cell 
lines in which GUI and/or GU2 were expressed whereas 
MVC5AC was not. 

Tissue expression patterns of MUC5AC and GU were also 
examined using a commercially available RNA panel of normal 
human tissues (Fig. 5B). Although several adult tissues such as 
stomach, lung, trachea, etc. obviously expressed MUC5AC 
mRNA, about half of the analyzed tissues were deficient in 
MUC5AC expression. On the other hand, GUI, GU2, and GU3 
were ubiquitously expressed in almost all the examined tissues. It 
should be noted that MUC5AC expression were not detected in 
adult brain, kidney, liver, placenta, skeletal muscle, spleen, 
thymus, thyroid gland, etc., despite the expression of the three 
GU homolog genes (Fig. 5B). 

Using non-malignant gastric tissues derived from the surgical 
specimens, we further evaluated the protein expression of 
MUC5AC and GUI by immunohistochemistry. As shown in 
Fig. 5C, Glil was weakly but ubiquitously expressed on gastric 
mucosa, which was diflFerent from the spotty immunostaining of 
MUG 5 AC. All the gastric tissues derived from the seven surgical 
specimens showed the same expression pattern. 

In total, expression patterns of MUC5AC and Gli clearly 
support our conclusion that GU is necessary for MUC5AC 
expression in human gastrointestinal cells. In addition, it is also 
obvious that expression of Gli is not enough for MUC5AC 
expression, suggesting that other mechanism should work coop- 
eratively on the transcription of MUC3AC gene. 

Multiple regulatory mechanisms work in the promoter 
region of MUC5AC gene 

To analyze the far distal promoter region oiMUC5AC gene, we 
generated reporter constructs containing 4010, 3000, 2000, and 
1433 bp upstream of ATG start codon and examined their 
promoter activities (Fig. 6A). Deletion from 4010 bp to 1433 bp 
region revealed gradual decrease in transcriptional activity, 
suggesting the presence of enhancer elements in this distal 
promoter region. However, even in the reporter construct 
inclucUng long 4010 bp, deletion of the 15 bp HCR-Gli region 
resulted in the drastic decrease of luciferase activity (Fig. 6A). 
Together, our results demonstrated that there are at least two 
regulatory elements in the MUC5AC promoter region: the 15 bp 
sequence between —111 bp to —125 bp which critical for the 
promoter activity, and the wider regulatory region in —4010 bp to 
— 1433 bp which can enhance the promoter activity. 

To further examine the possible epigenetic regulation on 
MUC3AC expression, four gastrointestinal cell lines and two 
glioblastoma cell lines were treated with demethylating agent 



(5Aza-dC) and/or HDAC inhibitor (TSA). Transcription of 
MUC5AC was slightly upregulated in all the four MUC5AC- 
deficient cell lines (Fig. 6B), suggesting that epigenetic regulation 
such as methylation and/ or histone acetylation additionally works 
as the regulatory mechanism oi MUC3AC expression. From the 
results of four MI/C5/4C-deficient ceU lines, methylation seems to 
play some universal role on suppression of MUC3AC expression. 
Synergistic effects of 5-Aza-dC and TSA in A172 and T98G cells 
also suggested that modification of histone acetylation may have 
some influence on MUC5AC transcription. 

Discussion 

MUG 5 AC, a secreted mucin highly detected in the superficial 
gastric epithelium, is a well-known gastric marker often used for 
classification of stomach cancer [20,39]. We have previously 
reported that expression of MUC5AC in stomach decreases in 
association with development of intestinal metaplasia [6] . In the 
present study, we further showed that MUC3AC expression is 
apparendy related to the tumor stage: advanced gastric cancers 
present reduced levels of MUG SAC compared with early gastric 
cancer. Despite the anticipated importance of MUC3A(^ regula- 
tion, the molecular mechanisms underlying MUC3AC expression 
in gastrointestinal cells remain poorly understood. In this study, we 
searched 5'-upstream of MUC3AC gene and identified the HCR- 
Gli at —125/— 111 bp in its promoter region. Our overexpres- 
sion/knockdown analyses and luciferase reporter assays revealed 
that Gli induced MUC3AC expression through the HCR-Gli in 
gastrointestinal cells. In addition, ChIP analysis showed that Glil 
directly binds to HCR-Gli. The results of RT-PCR using 30 
gastrointestinal tumor ceU lines and an RNA panel from systemic 
normal tissues suggested that GU is necessary but not sufficient for 
MUC3AC expression, which was consistent with the results from 
immunohistochemistry of clinical specimens. Furthermore, we 
found that the far distal upstream region from —1433 to — 
4010 bp enhances the Gli-dependent MUC3AC expression, at 
least partly based on the epigenetic modulation. 

Recendy, Inaguma et al. reported that Glil activates the 
MUC3AC promoter through two putative Gli-binding sites (GBSl 
and GBS2) in pancreatic cancer cells [25]. However, the 
significance of MUC3AC expression is different in gastrointestinal 
cells and pancreatic cells. Namely, abundant expression of 
MUC5AC is observed in normal epithelial cells of stomach but 
its expression is decreased accompanied with development of 
intestinal metaplasia [14,20]. On the other hand, MUC3AC is 
generally undetectable in normal pancreas tissue and is sometimes 
ectopically induced in pancreatic tumor cells. In fact, there are two 
dilferences between Inaguma's report and our results. First, 
Inaguma et al. showed that both GBSl and GBS2 respond to 
upregulation of MUC3AC expression by GU. The GBSl overlaps 
to our identified HCR-Gli, but our results indicated that the GBS2 
(corresponding to —159/ — 168 bp in Fig. 2) has no effect on 
MUC3AC expression (Fig. 3A). Second, Inaguma et al. reported 
that Glil expression correlates with MUG5AC induction by 
immunohistochemistry. However, our results revealed that 
expression of Gli does not always accompany with MUC3AC 
expression in gastrointestinal cells (Figure 5). These inconsistencies 
might be due to a dilference of originated tissues (pancreatic or 
gastrointestinal) . 
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Figure 5. Expression profiles of Gli and MUC5AC in various tumor cell lines and human normal tissues. (A, B) RT-PCR detecting MUC5AC 
and Gli homologs in the 30 tumor cell lines (A) and systemic normal tissues (B). (C) Immunostaining of IVIUC5AC (left panel) and Glil (right panel) in 
sequential sections of non-malignant gastric mucosa. 
doi:10.1371/journal.pone.0106106.g005 



No induction oi MUC5AC expression in AZ521 cells (Fig. 4A), 
which completely lack the SWI/SNF complex-based chromatin 
remodeling activity [28], may suggest the interaction between Gli 
and the SWI/SNF complex. Although there have been no 



previous reports concerning the interaction with Gli, SWI/SNF 
chromatin remodeling complex has been reported to bind with 
various transcription factors such as AP-1 [40], CREB [41], MyoD 
[42], Cdx [8], GATAl [43], and so on. 
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Figure 6. Multiple regulatory mechanisms of A^t/C5/ICgene expression. (A) Luciferase reporter analysis of a series of constructs covering the 
far upstream sequence in the MUC5AC promoter. Data represent the mean of luciferase activities in SW480, SH-10-TC, and KE-39 cells measured at 
24 h after transfection. The error bars showed the standard deviation of the results from three independent experiments. (B) Expression of MUC5AC 
analyzed by RT-PCR, using the 4 gastric cancer-derived and 2 glioblastoma-derived cell lines treated with 5-Aza-2'-deoxycytidine (5-Aza-dC) and/or 
trichostatin A (TSA). (C) Schematic representation of the presumed multiple regulations of human MUC5AC qene. 
doi:1 0.1 371/journal.pone.01 061 06.g006 



The 5'-upstream distal region from —1433 to —4010 bp in 
MUC3AC promoter showed obvious enhancing effect on 
MUC3AC expression (Fig. 6A). Interestingly, we noticed that the 
region from —1455 to —3187 bp contains 86 copies of tandem 
sequences 5'-TCA(C/T)TCA(C/T)-3' (Fig. 6C) [44]. The tandem 
repeats are conserved among mammals especially among primate 
and cetartiodactyla. Human, olive baboon, rhesusmonkey, Idll- 
erwhale, and dolphin have more than 60 repeats of the tandem 
repeats in their 4010 bp upstream of MUC5AC promoter region. 
In contrast, horse, dog, cat, panda, and hamster have only less 
than 10 repeats of the tandem repeats. We have no information for 
this tandem repeat sequences at present, but it might has some 
roles on regulation of MUC5AC expression. 

Recently, Yamada el al. showed that the CpG methylation 
status of around —3.7 kb in the MUC5AC promoter is associated 
with MUC5AC expression in some cancer cell lines, suggesting 
that the expression of MUC5AC is epigenetically regulated in the 
distal promoter region [45]. In fact, we found increase in 
MUC5AC expression by 5Aza-dC and/or TSA treatment. 



indicating some epigenetic mechanism involved in MUC5AC 
gene regulation (Fig. 6C). 

To understand gastric tumorigenesis mainly occurred in the 
atrophic mucosa with chronic gastritis and intestinal metaplasia, 
we beheve the disrupted balance between intestinal and gastric 
differentiation must be important [8,28]. For intestinal differen- 
tiation, many previous reports indicated that Cdx must be an 
indispensable key molecuk; based on the transcriptional regulation 
of many intestinal difiFerentiation marker genes [8,9,10,1 1,12]. On 
the contrary, the critical key molecule for gastric differentiation 
has not been identified yet. Despite the through screening, we have 
not found the essential regulator for our identified gastric marker 
Calhepsin E (data not shown) [6], and we have elucidated that 
typical gastric marker MUC3AC gene is partly regulated by Gli, a 
universal transcription factor which alone cannot realize the tissue- 
specific expression of MUC3AC gene. Unlike intestinal differen- 
tiation, we now speculate that gastric difiFerentiation may depend 
on cooperative regulatory mechanism of some universal transcrip- 
tion factors and epigenetic modulations. 
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In summary, the present study has shown that Gli regulates 
MUC5AC gene expression via direct protein-DNA interaction 
through the highly conserved 15 bp secjuence between —125 and 
— 111 bp in the promoter region of MUC5AC. Furthermore, 
immunohistochemical analysis and RT-PCR using systemic 
normal tissue revealed that Gli is necessar)^- but not sufficient for 
MUC5AC expression. We conclude that MUC5AC expression is 
regulated by combination of multiple regulatory mechanisms such 
as universal transcription factors and epigenetic modulations. 
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